ENERGY INEQUALITIES AND DISPERSIVE ESTIMATES FOR 
WAVE EQUATIONS WITH TIME-DEPENDENT COEFFICIENTS 
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Abstract. We consider wave models with lower order terms and recollect some re- 
cent results on energy and dispersive estimates for their solution based on symbolic 
type estimates for coefficients and partly stabilisation conditions. The exposition is 
complemented by a collection of open problems. 



Model problem. We consider on R< x R™ the Cauchy problem 

u tt - a 2 (t)Au + 2b(t)ut + m 2 (t)u = 0, (la) 
u(0, •) = ui, u t (Q, •) = u 2 (lb) 

for a wave equation with time-dependent coefficients and ask for the interplay between 
properties of coefficients (and their long-time behaviour) and large-time asymptotics 
of solutions (measured in terms of energies or certain norms). There is an extensive 
literature about results; we will sketch some of them and focus on related open problems. 

Problem (1) looks far more general than it is. Using Liouville transforms it is possible 
to reduce the problem to equations with constant a(t) and vanishing b(t). We will not 
pursue this further and merely restrict our consideration to equations with variable speed 
and equations with constant speed and variable dissipation. 

I. Some classical results on bounded variable propagation speed. 

I. (a) For the sake of completeness we start our consideration with the wave equation 

u u - An = 0. (2) 
Then it is well-known that the corresponding energy 

EM) = X - J (|Vu(t,x)| 2 + \u t (t,x)\ 2 )dx (3) 

is preserved, K(u;t) = E(u;0). Furthermore, if we replace the L 2 -norm by an L 9 -norm, 
q > 2, we gain decay. This is known as Strichartz type decay estimate or dispersive 
estimate and goes back to classical papers of Strichartz, [Str70], Pecher, [Pec76], and von 
Wahl, [vW71]. It plays an essential role in considerations of nonlinear wave equations. 
In a precise formulation the result reads as 

||Vu(t s -)IU« + < ^(l + ^-^^-i^lKDp+Vllw + \\(B)^u 2 \\ LP ) (4) 

for p£ (1, 2], pq = p + q and r p = n(l/p — 1/q). As usual we denote (D) = (I — A) 1 / 2 . 
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There exist several alternative formulations of this estimate, we want to point out two 
of them. First, if we replace (1 + t) by t (and thus allow the decay function to be- 
come singular in t = 0) the regularity can be sharpened to r p = (n + l)(l/p — l/q)/2. 
Furthermore, this statement is then equivalent to the Strichartz' estimate 

\H\LjLr + ||Vll|| £ oo £ a + ||u t || L oo L 2 < C(||Vni|| L 2 + \\U 2 \\ L 2) (5) 

with 2 < q < oo, 2 < r < oo, l/q + n/r = n/2 and 1/q + (n - l)/(2r) < (n - l)/4 as 
pointed out in [KT98]. 

The simplest case of time-dependent propagation speed is a periodic coefficient, 

u tt -a 2 {t)Au = 0, a(t + T) = a(t) (6) 

with a E C 1 (IR), < a\ < a{t) < a 2 and dta ^ 0. Regardless of the 'size' of the 
coefficient, this destroys any reasonable energy estimate as pointed out by Yagdjian, 
[YagOl]. We will give a short outline of the argument. Applying a partial Fourier 
transform to (6) yields Hill's equation 

u tt + Xa 2 (t)u = (7) 

with spectral parameter A = |£| 2 £ [0, oo). Due to a variant of Borg's theorem (see, 
e.g., the text book of Magnus- Winkler, [MW66] or Colombini-Spagnolo [CS84] for a 
refinement) there exists an open interval X C [0, oo) such that for all A € X equation (7) 
has an exponentially increasing solution. Taking appropriate Cauchy data for (6) with 
Fourier support on X we can construct solutions with exponentially growing energy. 

Theorem 1 (Yagdjian, [YagOl]). Assume a G L ; 1 oc (lR) is non-constant, periodic and 
a.e. positive. Then there exist Cauchy data u(0, •) = u\ and u t (0, •) = u 2 from Schwartz 
class such that the energy of the solution to (6) satisfies 

^logEM)^ 

t^oo logt 

One can say a bit more about the structure of solutions to (7). There exist nested 
sequences A^, = A^ < AT/ < < < A^ < • • • , such that for all A € (A+, A~ +1 ) 
solutions are quasi-periodic, while for A G Xj = (\J,\f) solutions are combinations 
of exponentially increasing and exponentially decreasing contributions. Some of the 
intervals Ij may be empty, but not all. 

1.(7) We can not ask for conservation of energy if we consider variable propagation speed 
(and no compensating lower order terms). However, we may ask for conditions on the 
coefficient function a(t) such that the solutions to 

u tt -a 2 (t)Au = 0, u(0,-)=«i, u t (0,-)=u 2 (9) 

satisfy Strichartz type decay estimates (and therefore also uniform energy bounds). A 
first answer was given by Reissig-Smith in [RS05]. We assume 

< ai < a{t) < a 2 , (10a) 
d k ( 1 \ k 

^a(t) <C fc f— J , fc=l,2,...,n+l. (10b) 
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Then the following statement holds true: 

Theorem 2 (Reissig-Smith, [RS05]). Assume a G C n+1 (IR) satisfies (10). Then the 
solution to (9) satisfies the a priori estimate 

||Vn(t,-)||^ + lkt(*,-)ll^ < C^Cl+^-^^-^dKD^+^xtl^ + IKD^^H^) (11) 
for p G (1,2], pq = p + q and r p = n(l/p — 1/q). 

If we are interested in the energy alone, the constructed representation of solutions im- 
plies the two-sided energy estimate C\E(u; 0) < E(u; t) < C2E(u; 0) for certain constants 
C\ and C2. For this result, (10b) suffices for k = 1, 2. 

Condition (10b) can be weakened by the introduction of further log-terms (but increasing 
also the necessary number of derivatives). Then the estimate (11) loses a small, resp. 
finite, amount of decay. For details we refer to [RS05]. Results are sharp, there exist 
counter-examples to the estimate (11) if the conditions are slightly violated. 

Example 1. We give an example for an admissible coefficient function. The theory of 
Reissig-Smith allows to consider a(t) = 2 + sin(logi). For this coefficient the solutions 
to (9) satisfy Strichartz-type decay estimates. 

The proof of Theorem 2 is based on an explicit construction of the representation of 
solutions in terms of the coefficient function a(t). Properties of the representation are 
different in different parts of the phase space x M.^, within the hyperbolic zone {(1 + 
S> 1} it can be shown that solutions are represented by Fourier integrals 

VW) J 6^*1*1 ^WA(*,0/(0d£, (12) 

/ given in terms of the initial data and A being a symbol of order zero uniform in t. 
The construction is based on a diagonalisation argument within the hyperbolic zone. 

1.(6) The 'gap' between (0) and (7) can be partially filled. Hirosawa, [Hir07], weakened 
the symbol-like condition (10b) to smaller improvement per derivative in combination 
with a new assumption: the stabilisation condition. We consider the Cauchy problem 
(9) and assume for the coefficient 

< ai < a(t) < a 2 , (13a) 
rt 

\a(s)-a QO \ds<C(l + t) q , (13b) 





d k 

-—rQ,(t) 

dt k v ' 



^ Ck {—t) ' k = l,2,...,m. (13c) 
with m > 2, q <G [0, 1) and p > q + (1 — q)/m. Then the following statement holds true: 

Theorem 3 (Hirosawa, [Hir07]). Assume a G C m (M.) satisfies (13). Then there exist 
constants C\, C2 > such that the generalised energy inequality 

CiE(u; 0) < E(u; t) < C 2 E(n; 0) (14) 

is valid. 
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On the contrary, there exists a coefficient satisfying (13) with p < q for which no uniform 
energy bound holds true. Therefore, also this result is sharp. 

Example 2. a(t) = 2 + sin(i a ), a > 0, does not satisfy (13). The stabilisation condition 
is violated. Indeed, with doo = 2 we obtain J* * | sin(s Q )|ds = ^ J * | sm(6)\6 1 / a ~ 1 d9 ~ t. 

Example 3. To obtain an admissible coefficient we follow [HW09] and choose sequences 
r)j, Sj and tj, tj + Sj < ij+i, together with a function ip G Cq°(M) subject to supp ip C 
[0,1], -1 < 4>(t) < 1 and Jj |^(t)|dt = 1/2. Then 

^) = l + f>V>(^) (15) 

satisfies (13) provided that rjj < 1, Yl^iVj^j < ^j+i an< ^ *j — Such sequences 

exist, e.g., tj = 2^', Sj = 2?« and rjj = 

The proof of Theorem 3 is again based on explicit constructions of solutions. There 
are two major differences to the situation of Thm. 2. The hyperbolic zone is smaller 
{(l + t) 9 |£| 3> 1} (which makes it necessary to invoke a new argument for small fre- 
quencies) and the weaker assumptions for derivatives make it necessary to perform more 
diagonalisation steps. In consequence, representations by Fourier integrals contain an 
inhomogeneous phase function. 

Open Problem 1. Assume (13) for p > q and all k. Is is possible to derive Strichartz 
type decay estimates under these assumptions? 



II. Equations with increasing speed. Problems with increasing speed have been 
considered by Reissig-Yagdjian, [RY00]. Following their approach we write the coefficient 
function a(t) as 

a(t) = \(t)u(t) (16) 
with a monotonously increasing function A(i) and a (bounded) oscillating part uj(t). We 
denote a primitive of \(t) as A(i) = 1 + J Q * A(s)ds and assume a two-sided estimate 

\'(t)A(t) ~ (X(t))\ (17) 

We refer to such functions A(i) as admissible shape functions and investigate the Cauchy 
problem 

u tt - X 2 (t)u 2 (t)Au = 0, u(0,-) = u u u t (0,-)=u 2 (18) 
for suitable u(t). 

II. (a) We shortly recall the result of Reissig-Yagdjian, [RY00], again neglecting loga- 
rithmic terms and the related classification of oscillations to simplify statements. We fix 
an admissible shape function A(i) and assume 

< ciA(t) < a(t) < c 2 X(t) (19a) 

^a(t)| <C k \(t)(j^ , fc = l,2,...,n + l, (19b) 
Then the following Strichartz type decay estimate holds true: 
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Theorem 4 (Reissig-Yagdjian, [RYOO]). Assume a G C n+1 (R) satisfies (19). Then the 
solution to (18) satisfies the a priori estimate 

||A(t)Vu(t,-)|| L , + K(t,-)|| L , <^yAW(A(t))- 2 ^ i( i-| ) (||(Dp+V||L P + IKDr^ 2 || L p) 

(20) 

forp G (1, 2], pq = p + q and r p = n(l/p — 1/q). 

If we consider the energy estimate alone, this yields E\(u;t) < A(i)(||iti||^i + ||u2||^ 2 ) 
for the adapted energy 



E x (u; 



t ) = \j (|A(t)Vn(t,x)| 2 + \u t (t,x)\ 2 )dx. (21) 



II.(/3) Again, we want to weaken the conditions to allow stronger oscillations in combina- 
tion with a stabilisation condition in the spirit of [Hir07]. We assume that a(t) = X(t)u(t) 
satisfies 

= , A » (() |< A(t) (A|) 2 , li=p ~ <2 , (22a) 

< ci < w(t) < c 2 , (22b) 
the stabilisation condition 

f X(s)\uj(s) - cJoo|ds < C0(t) = o(A(t)) (22c) 

for an auxiliary function 0(t) < A(i), 0(0) = 1, and a constant cJqo, the symbol- like 
estimate 



-rra(t) 
dt k y 1 



<C k X(t)(E(t))- k , fc = l,2,...,m (22d) 

with an (increasing) function S(t), A(i)E(t) > C0(t), and 

(A(s)) 1 ~ m (~(s))-"Ms < C(0(t)) 1 - m . (22e) 
Under these assumptions the following energy estimate holds true: 



Theorem 5 (Hirosawa-Wirth, [HW09]). Assume a G C m (R) satisfies (22). Then solu- 
tions to the Cauchy problem (18) satisfy the energy estimate 

E A (u;t) < \(t)(\\ Ul \\ H i + \\u 2 \\v). (23) 
Furthermore, if (ui,U2) ^ (0,0) the limit 

lim -i-E A («;t)^0 (24) 
t^oo A(rJ 

exists and is non-zero. 

Open Problem 2. Derive dispersive estimates under similar assumptions. This should be 
a consequence of the solution to Problem 1. 
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III. Equations with non-effective dissipation. We consider the influence of dissi- 
pation terms and investigate 

u a - Au + 2b(t)u t = 0, u(0,-)=«i, u t (0,-)=u 2 . (25) 

This Cauchy problem is related to (18) with increasing propagation speed by a change 
of variables. For completeness we give the correspondence. Let v(t,x) be a solution to 
vtt — \ 2 {t)Av = 0. Then u(t, x) = v(A~ 1 (t),x) — with A _1 (i) the inverse function to 
A(t) — satisfies 

vt(&-Ht),x) 

Ut= A(A~ 1 (t)) ■ (26a) 

W « = ^(A-ift)) " *MoP (A ®> X) = AU ~ (26b) 
which is of the form (25); X(t) = (1 + t) e corresponds asymptotically to 2b(t) ~ j^ijifi- 

III. (a) In [Wir04], [Wir05], [Wir06] the author gave an overview on energy and disper- 
sive estimates for solutions to (25) under assumptions related to [RY00], [RS05]. We will 
quote a corresponding result. We assume for the coefficient function 

b(t) > 0, (27a) 
lim sup io(i) < ^ (27b) 

t^oo 2 
d k / 1 \l+k 

^b(t) <C k l — ) , fc = l,2,...,n + 2. (27c) 
Then the following Strichartz type decay estimate holds true: 

Theorem 6 (Wirth, [Wir06]). Assume b £ C n+2 (R) satisfies (27). Then the solution 
to (25) satisfies the a priori estimate 

iiv«(t, oiiw + k(*, oiu- < ^^(i +r V( ^ ) (ii(D> rp+ viiL, + ii(D)^u 2 || LP ) 

(28) 

/or p G (1, 2], = p + q and r p = n(l/p — 1/q) and with the auxiliary function 

f3{t) = exp f b(s)ds. (29) 
Jo 

Furthermore, for non-zero initial data {u\,U2) € H 1 x L 2 the limit 

lim fi 2 (t)E(u;t) + (30) 

t— >oo 

exists and is non-zero. 

In particular, (30) can only be true, if tb(t) is asymptotically small. As pointed out by 
Hirosawa-Nakazawa, [HN03], the example b(t) = +*) with fi > 1 satisfies 

lim t 2 E(u; t) = (31) 
for all data and, furthermore, the exponent 2 is best possible (see, e.g., [Wir04]). 
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III. (/?) The conditions (27) can be weakened in the spirit of [Hir07] to allow faster 
oscillations of the coefficient. We assume that 2b(t) = fi(t) + a(t), where fi(t) is a 
monotonically decreasing shape function and a(t) carries the oscillations. Conditions 
are similar to II. (/?) and follow [HW08]. We assume 



fj,(t)>0, fi'(t)<0, limsupi/i(t) < 1, 



sup 

t>0 



/ cr(s)ds < oo, 
Jo 



cxp 



<r(s)ds 



with an auxiliary function @(t) = o(t), t — > oo, 
d fc 



< C fc (E(t)) 



-fc-1 



d# < ce(t), 



fc = 0,1,2, 



, m 



(32a) 
(32b) 
(32c) 

(32d) 



with E(t) > 0(t) and 

(£(*))— Ms < (9(t))- m . 

Then the following statement holds true: 

Theorem 7 (Hirosawa-Wirth,[HW08]). Assume b G C 
energy of the solution to (25) satisfies 

E(u;t) < C^Tj (INIIff 1 + IWIi^ 



(32e) 

satisfies (32). TTien i/te 
(33) 



wi/j /3(i) = exp J** 6(s)ds ~ exp f*[i(s)ds. Furthermore, for non-zero initial data the 
limit 

lim /? 2 (t)E(u;t) / 



t^oo 



(34) 



exists and is non-zero. 



Example 4. We can use fi(t) = and a(t) = fi(t) sin(t Q ) with a E (0,1). Then all 
conditions are satisfied. Note, that the corresponding approach due to Reissig-Smith 
allows only sin(logt) terms. 

Example 5. We can come closer to sin(t) by using fi(t) = k>g(e+t) an< ^ a (f) = 
/i(t)siii(t/log(e + t)). 

Note, that a(t) has no influence on the decay rates, only on appearing constants. a(t) 
can not be seen as a small perturbation of n(t); it might be arbitrary large. 



IV. Equations with effective dissipation. Following the classification of [Wir05], 
[Wir07], we speak of effective dissipation if the lower order term 2b{t)ut changes the 
large time asymptotics of solutions to (25) in an essential way. 
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IV. (a) The damped wave or telegraph equation 

%-A« + u ( = 0, u(0,-)=ui, u t (0,-)=u 2 (35) 

was studied by Matsumura in [Mat 76]. He investigated semilinear perturbations of this 
equation based on 

Theorem 8 (Matsumura, [Mat76]). The solutions to the damped wave equation (35) 
satisfy the a-priori estimates 

\\d^u(t,x)\\ L 2 < C(l + t)~^- k -^ (\\ Ul \\ Hk+la i + ||« 2 ||h*+M-i 

+ + ||u 2 || L i) (36a) 

\\djfd%u(t,x)\\ L °° < C(l + t)~ l ^~ k ~^ (||ltl|| H r+fc+|a| + ||U2||ffr+fc+|a|-l 

+ ||ni|| L i + ||u 2 ||li) (36b) 

with r > . 

The estimates are different in structure compared to the Strichartz type estimates for the 
free wave equation and the related estimates discussed so far. While the used regularity 
is related to Sobolev embedding (and therefore hyperbolic in nature) the decay rates 
correspond to estimates for the heat equation (as, e.g., used by Ponce in [Pon85]). 
This parabolic structure becomes even more apparent in the diffusion phenomenon due 
to Nishihara, [Nis97], [Nis03], and Yang-Milani, [YMOO] or the estimates of Narazaki, 
[Nar04]. 

Theorem 9 (Nishihara, [Nis03]). Let u(t,x) be a solution to (35) in R t x M%., w(t,x) 
the corresponding solution to 

w t = Aw, w(0, •) = u\ + u 2 (37) 

and v(t,x) a (suitable) free wave. Then the L p -L q estimate 

\\u(t, •) - w(t, •) - e-3v(t, < C(l + ty^p-v ) ~ 1 (||«i||lp + ||n 2 ||Lp) (38) 
holds true for all 1 < p < q < oo. 

The decay rate is by one order better than expected for both \\u(t, -) and \\w(t, -)IU 9 
(while the regularity of the data on the right hand side is not sufficient to obtain these 
decay estimates — singularities are eliminated by the free wave v(t,x)). 

For both statements it is essential that small frequencies determine the asymptotic be- 
haviour of solutions. If we localise to high frequencies (or to a hyperbolic part of the 
phase space), we obtain exponential decay rather than the polynomial decay given in 
the theorems. 

IV. ((3) In [Wir07] the author discussed variable coefficient versions of these results; we 
quote a slightly simplified version (c.f. [RW06] or [Wir05] for this formulation). Let b(t) 



ENERGY INEQUALITIES 

be subject to 



b(t) > 0, b'(t) ^ 0, lim tb(t) = oo, (39a) 

t— +oo 



d k 



<C fc 6(t)( — ) , fc = l,2. (39b) 



1 
1 

Then the following decay estimate holds true: 

Theorem 10 (Wirth, [Wir05]). Assume b G (7 2 (E) satisfies (39). TTten i/te solution to 
(25) satisfies the a priori estimate 

Ht,-)\\ Lq <C pqr il + J ^yj (IKD^uxUlp + IKD^-^Ulp) (40a) 

||Vu(t,-)|| LB <C pqr (l + J^ 2 (||(D)^+V|| iP + ||(DPn 2 || iP ) (40b) 

K(V)I|l* <C pqr (l + J -j^j 2P " (IKD^+Vllix + IKD^uallip) (40c) 
for l<p<2<g<oo and r p > n(l/p — 1/q). 

The proof of this statement is constructive in the sense that the main terms in the 
representation of solutions are given explicitly. The assumption tb(t) — > oo guarantees 
that the high frequency part decays faster than the low frequency part and is essential for 
the parabolic type estimates. Under slightly more restrictive assumptions the diffusion 
phenomenon can be obtained. We refer to [Wir07] for the details. 

IV. (7) We conclude the overview with some remarks concerning periodic dissipation 
terms. In contrast to periodic propagation speed this does not destroy Matsumura's 
decay estimates nor the diffusion phenomenon. The following two statements are based 
on [Wir08]. 

Theorem 11 (Wirth, [Wir08]). Assume b € ACi oc (R) is absolutely continuous, positive 
a. e. and periodic. Then the solution to (25) satisfies the a priori estimate 

\\u(t,-)\\ Lq <C pqr (l + t)-^--*\\\(-Dy*> Ul \\ LP + ||(D)^- 1 u 2 || L p) (41a) 

||Vu(t,-)|| LB <C; ?r .(l + t)-t^-i)-i(||(Dp + V||L, + ||(Dp« 2 || LP ) (41b) 

||«t(*,-)IU« < q^(l + t)" 5( '~* )_1 (ll<D> rp+1 «itl^ll<D> r ^lli>) (41c) 
for l<p<2<q<oo and r p = n(l/p — 1/q). 

Open Problem 3. Close the gap between IV. (f3) and IV. (7), i.e., develop energy estimates 
based on less regularity (and/or less improvement per derivative). 

Theorem 12 (Wirth, [Wir08]). Assume b G ACi oc (R) is absolutely continuous, positive 
a. e. and periodic. Then there exist constants a and (3 depending only on the coefficient 
function b such that the solution to (25) and the corresponding solution to 

wt = aAw, w(0, ■) = u\ + (5u2 (42) 
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satisfies 

\\u(t, •) - w(t, -)\\ L 2 < C(l + ^^(llmll^i + ||n 2 || L2 ). (43) 

The constants a and (3 are explicit and given in [Wir08]. A similar estimate holds true if 
we replace the L 2 -norm on the right hand side by an L 9 -norm, q > 2, or measure spatial 
derivatives of the difference. 

Both theorems are essentially based on the idea that small frequencies determine the 
asymptotic behaviour. Although the periodic coefficient leads to a variant of Hill's 
equation after partial Fourier transform and a change of variables, we only need to 
understand the neighbourhood of £ = 0. It can be shown that away from this exceptional 
frequency, solutions decay exponentially, while nearby (which is inside an instability 
interval) the structure of the fundamental solution is known explicitly due to Floquet's 
theory (cf., [MW66] or [Wat22, Chapter IXX]). 

V. Concluding remarks. We focussed on the influence of variable speed and variable 
dissipation. There are further results on the treatment of variable mass terms, see e.g. 
the treatment of Bohme- Reissig, [BR09] , of scale invariant mass or [HR03] on interactions 
between variable mass and variable speed. 

There is a fundamental difference between oscillations in the principle part and oscilla- 
tions in lower order terms; mainly due to asymptotic relations to problems of different 
type. While in the first case sharp results are available (within certain symbol-like 
classes of coefficients), the example of effective dissipation shows an essential gap of 
understanding between [Wir07] and [Wir08], which needs to be closed in the future. 
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